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Abstract: The Minimal Supersymmetric Standard Model (MSSM) can include two 
soft breaking terms which are often neglected: a non-analytic scalar trilinear coupling 
and a Higgsino bilinear term. A set of high-scale boundary conditions consistent 
with the reparameterisation invariance which the model possesses is obtained. The 
three-family renormalisation group equations for the MSSM with these terms are 
presented. The ranges of the universal high-scale values of these couplings which lead 
to an acceptable TeV-scale theory are obtained, as is the supersymmetric particle 
spectrum at this scale. The effect of the new terms on fine-tuning is presented. 
SOFTSUSY, an existing program for calculating SUSY particle spectra, has been 
used, with as few modifications as possible. 
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1. Introduction 

A possibility for new physics beyond the standard model is supersymmetry (SUSY) [1, 
2]. If fermionic generators are added to the bosonic generators of the Lorentz group, 
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the new space-time symmetry is supersymmetry. In a theory with exact supersym- 
metry, all particles have a partner of equal mass but opposite spin-statistics. Can- 
cellations between bosonic and fermionic loops prevent radiative corrections from 
driving scalar masses up to the highest scale present, assumed to be the GUT or 
Planck scale, 10 16 to 10 19 GeV, solving the naturalness problem of the Standard 
Model. In addition, the renormalised electromagnetic, weak, and strong couplings 
can be made to converge to an approximately common value at the grand unification 
scale. 

Since supersymmetry is not observed amongst the already discovered particles 
it must be a broken symmetry. The superpartners of the observed particles must 
be undiscovered particles. The simplest possible SUSY extension of the Standard 
Model, with a new superpartner field for each standard model field, and the addition 
of a second Higgs scalar doublet, is called the minimal supersymmetric Standard 
Model, or MSSM. 

In the MSSM, the supersymmetry breaking terms are often assumed to be 
flavour-independent and/or gauge-factor independent at the high scale, and to split 
as they evolve to low scales under the renormalisation group equations (RGEs). The 
model with this assumption is called the constrained MSSM (CMSSM) or minimal 
supergravity (mSUGRA). The universal nature of the SUSY breaking terms is mo- 
tivated by the idea that supersymmetry breaking is mediated by some flavour-blind 
particle such as the graviton. 

Within mSUGRA, the values of fundamental SUSY-breaking parameters are im- 
posed as boundary conditions on the running SUSY-breaking masses and couplings 
at high scale, usually taken to be M x = 10 16 GeV. Experimental masses of su- 
perparticles are obtained by evolving the Lagrangian to the weak scale using the 
renormalisation group equations. A solution consistent with the mSUGRA bound- 
ary conditions at Mx and the known Standard Model constraints at Mz must be 
found. This constitutes a boundary value differential equation problem, and can be 
solved numerically. 

The MSSM Lagrangian is usually claimed to include all possible "soft supersym- 
metry breaking" terms, terms which split the masses and couplings of particles and 
their superpartners, but which do not remove the supersymmetric protection against 
large radiative corrections to scalar masses. It is also supposed to include all possible 
SUSY conserving terms given the particle content. 

In fact, the "R-parity violating terms", which, if unconstrained, allow proton 
decay and flavour changing neutral currents, are also possible, but conventionally 
excluded. If these terms are excluded, superparticles are created and destroyed only 
in pairs, giving rise to the usual SUSY phenomenology, where the lightest supersym- 
metric particle is stable. 

There exists a further set of possible additional "non-standard" soft supersymme- 
try breaking (NSSB) terms allowed by the symmetries, as remarked upon in [1,3-7]: 
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a non-analytic scalar trilinear coupling, and a Higgsino bilinear coupling: 

c nssb = c^H^QiUj + CfjHZQidj + CfjHZL&j + jiH x H 2 + c.c, (1.1) 

where i,j are family indices, and weak isospin, spinor, and colour indices are sup- 
pressed. C{j and jl are the new couplings. 

In [3] the C terms are written down and then taken as zero. In [4] Feynman 
rules and mass matrices for these terms are presented. The non-standard terms 
are written down in [1], and including them in SUSY studies is advocated there. 
In [5] and [6] the NSSB terms are used to generate flavour mass hierarchies through 
radiative corrections. 

In [7] 1-loop /3-functions are obtained for a SUSY model with NSSB terms, 
and a general particle content and gauge group. These are then specialised to a 
one-family MSSM. The region of m ,tan/3 parameter space where an acceptable 
electroweak vacuum results is obtained, for a given value of fl(Mx) and with /x = 0. 
Co = C u /Y u = C d /Y d = C e /Y e is allowed to vary in order to obtain an acceptable 
vacuum, as discussed in section | of this paper. When Co = ft the theory is super- 
symmetric (C and jx can be removed by a reparameterisation as discussed in section ^ 
of this paper.) This is a fixed point of the RGEs, and its stability is discussed in [7]. 
In [8] those authors present 2-loop NSSB /5-functions for general particle content and 
gauge group, and further discuss the fixed-point properties of the NSSB RGEs. 

In [9] the impact of NSSB terms on flavour- changing neutral currents in the 
MSSM is considered. In [10] the NSSB terms are included together with dimension- 
less higher-order scalar couplings to examine the origin of intermediate scale physics. 
In [11,12] the NSSB terms are used in models of neutrino masses and rare processes 
in top decay. In [13] the NSSB terms are used to fix the correct mass of the 6-quark 
in a model where the decay b — > 57 is generated radiatively. 

These previous works have not considered how the mSUGRA assumptions can 
be adapted to be used with NSSB, determined the size of the valid NSSB parameter 
space for fixed values of the standard MSSM parameters, presented mass-spectra 
for the NSSB MSSM, or calculated the fine-tuning of the model. These issues are 
examined in the present work. 

In the following sections the MSSM with NSSB is considered in detail. Section |2] 
is a review of the reasons why the NSSB terms have previously been neglected. In 
section numerical methods which can be used to simulate the renormalisation group 
equations are discussed. (For more detail on this question, consult appendix [B].) In 
section |4] a reparameterisation invariance is introduced, which the inclusion of these 
terms gives to the MSSM. In section [5] we develop a version of the "constrained 
MSSM" or "minimal supergravity" (mSUGRA) high-scale boundary conditions con- 
sistent with the reparameterisation symmetry. In section |] the calculation of the 
3-family renormalisation group equations (RGEs) for this model is considered. These 
are presented in appendix |A]. The impact of these terms upon electroweak symmetry 
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breaking is discussed in section [7]. In section ||| an interesting subset of the particle 
masses for a standard mSUGRA point but with non-standard breakings is presented. 
In section R] the reasons for the limits on the magnitudes of the non-standard cou- 



plings are explored. In section 10 numerical measures of the fine-tuning of the model 
are considered, and results for the fine-tuning presented and discussed. Appendix [A| 
lists the NSSB contribution to the 1-loop renormalisation group equations for the 
model. Appendix [B] examines the numerical methods used in this work. 



2. Why have the NSSB terms traditionally been excluded? 

The terms of eq. have traditionally been neglected because they are thought 

to lead to large radiative corrections when there are scalar fields present which are 
uncharged under the entire gauge group. 

This argument cannot be used in the context of the MSSM, which contains no 
singlets at the SUSY scale. In a model with some higher mass singlet field of mass 
m s , together with a very massive field of mass mx, the correction, (see figure |l|) to 
the relation between the Higgs pole m^ p and running m^ r masses is: 
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where C s is the coupling between the Higgs and the singlet, 
and Cx is the coupling between the singlet and the heavy 
field. S4 is a collection of dimensionless constants. If the 
singlet field has a mass significantly below the high scale, 
the theory is thus unstable. However, if any singlet fields 
have masses of order Mx, such as in SUSY desert scenarios, 
then there is no problem with including NSSB terms. 

3. Numerical calculation of SUSY particle 
spectra 



m. 



m. 



Figure 1: The scalar 
singlet tadpole correc- 
tion to the Higgs prop- 
agator. 



In SUSY models, the supersymmetric Lagrangian terms are fixed at the experimental 
scale Mz from Standard Model particle data. The supersymmetry breaking terms 
are fixed at some high scale, Mx, according to a scheme such as mSUGRA. 

The high- and low- scale boundary constraints, together with the renormalisation 
group equations, and the intermediate scale radiative electroweak symmetry break- 
ing (REWSB) constraint (see section ^), form a system of differential equations; a 
"boundary value problem" . 

The numerical methods used to implement this procedure are discussed in detail 
in appendix O. Three such methods are discussed. The first is the formal "shooting 
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method" where the constraint equation at one boundary is solved for those param- 
eters which are unconstrained at the other boundary. This method cannot be used 
for SUSY, because of the complex form of the boundary conditions at the low scale. 

The second method, which may be called the "drift method" involves alternate 
imposition of each set of boundary conditions. This constitutes a recurrence relation 
for successive values of those parameters unconstrained at one of the boundaries. 
The solution consistent with both sets of boundary conditions is a fixed point of the 
recurrence relation, but this point is not guaranteed to be stable or unique. If we 
assume the point is both stable and unique, however, simple repeated imposition 
of the boundary conditions, alternately at Mz and Mx, will result in convergence 
toward the required solution. This is the method used in ISASUGRA [14], and in 
SOFTSUSY [15], the code which was modified for this work. The C++ inheritance 
structure of SOFTSUSY permitted large parts of the code of SOFTSUSY to be 
re-used. 

Because the drift method becomes unstable for certain values of the NSSB pa- 
rameters, a third approach has been developed, a "single-ended" shooting method, 
where the fixed point of the recurrence relation is solved for, using a Newton-Raphson 
method. This is further discussed in appendix B. 



4. The NSSB Lagrangian and ^-reparameterisation 

The nonstandard SUSY breaking terms to be added to the MSSM Lagrangian are 
given in eq. (|1.1|) . Using the usual formula [16] for extracting the full Lagrangian 
from the superpotential, the following terms are those involving m 2 H , the Higgs 

1/2 

soft masses, /x, the Higgs bilinear superpotential term, and the NSSB couplings: 
C = ... +(/2 - /i)#i# 2 + (/i 2 + m^li/il 2 + (/i 2 + m 2 H2 )\H 2 \ 2 

+c.c, (4.1) 

where e are the Yukawa couplings. 

Therefore, the Lagrangian is invariant under the following reparameterisation: 

y, — > n + 5; 
fj, — > jj, + 5; 
m 2 H — ► m 2 H — 2y5 — 5 2 ; 

C[, ■ YjS. (4.2) 

This was partially pointed out in [7] where it was noted that one could write an 
MSSM Lagrangian (p, = 0, CL = 0, \i ^ 0) in terms of an NSSB model with zero /i 
but non-zero C(a and jl. Here a different emphasis is used:, the Lagrangian contains 
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all terms allowed by the symmetries, and the reparameterisation is used like a gauge 
freedom, to allow calculations to be made in the most convenient "gauge" . 

The above Lagrangian contains all terms involving ll, provided the Higgs bilinear 
breaking term is written as 

C = ...+7T%H 1 H 2 (4.3) 

and not 

£ = ...- BliH x H 2 (4.4) 

as is sometimes used, factoring out the ll. In the latter case, B would become a 
/x-reparameterisation dependent quantity. In order that this work can easily make 
contact with MSSM work conducted using B, we instead have the freedom to define 
B using the Lagrangian term: 

C=...-B(n-fi)H 1 H 2 . (4.5) 

On the MSSM subspace of parameter space, where jx and C{j are zero, this definition 
is consistent with eq. (|4.4j ), but B remains /i-reparameterisation invariant. 

5. The high-scale boundary conditions 

As mentioned in section |3| the supersymmetric Lagrangian terms are fixed at the 
experimental scale Mz from Standard Model particle data. The supersymmetry 
breaking terms are fixed at some high scale, Mx, according to a scheme such as 
mSUGRA. 

The exceptions to this are the superpotential /i-term, and the SUSY breaking 
term B. The magnitude of the Higgs vacuum expectation value (VEV) v 2 = vf + v\ , 
is known from the Z boson mass. In addition, a choice for tan/? at Mz is usually 
made. The values of ll and B fix both v and Vi/v 2 = tan/3, the ratio of the two 
Higgs vacuum expectation values, by minimisation of the Higgs effective potential 
(see section [7]) which is a known function of the model parameters. Thus, the values 
of ll and B cannot be imposed as a boundary condition, but are allowed to vary to 
fit the Higgs VEVs. 

While the above approach is the conventional one for MSSM calculations, in [7] 
the authors chose to use the reparameterisation freedom to take ll = 0, solving the 
SUSY /i-problem. Instead of leaving ll and B at Mx unconstrained, so that they can 
vary to fit Mz and tan/3, in this approach C and B are left unconstrained at Mx- 
(With C e , C d , C u taken as equal at Mx and jl at Mx given as an extra boundary 
condition.) 

In this investigation, we want to impose values for both C and jl as boundary 
conditions at My, with the same kind of universality assumptions as in mSUGRA, 
i.e. Oh = Kta. The pa_ „ w UI remain in the m0 del, and win, tog ethe r with 
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B, be varied to fit tan (3 and v. Since we do not choose fj, = the model retains its 
reparameterisation symmetry. 

This approach is perhaps preferable because it implements non-standard super- 
symmetry breaking with minimal changes to mSUGRA. However, there is a subtle 
difficulty, because we wish to use the mSUGRA assumption of universal high-scale 
soft scalar masses: the mass term for each of the sleptons, squarks and Higgs bosons 
should be the same at Mx- This assumption is inconsistent with the reparameter- 
isation invariance. A reparameterisation changes the Higgs soft masses but not the 
squark and slepton masses, breaking the equality To make the boundary condition 
well defined, we must therefore select the reparameterisation in which we want the 
scalars to be universal. 

As will be seen in section |6], the renormalisation group running will take place 
in fj, = reparameterisation. If we simply impose the mSUGRA boundary condi- 
tions in that reparameterisation, including a value for jl, we cannot have \x as an 
unconstrained parameter at high scale. 

Consider a model with m , Mi/ 2 , Co and jl at Mx given as constraints, and 
values for /i (us) an d B (Bg) and for the other parameters unconstrained at high 
scale, but which are consistent with the low energy constraints. (This discussion is 
illustrated in table [l|). On changing to the fi = reparameterisation used to run the 
RGEs, or the jl = reparameterisation used for some parts of the calculation, the 
values of \i and jl are mixed up, and only (jl — ji) is known. In addition, the Higgs 
soft mass-squareds become different from those of the squarks and sleptons. 

If we then run the RGEs, do the rest of one cycle of the calculation, and then 
return to Mx, we do not know how to separate fj, and ji. It seems we cannot return 
to the reparameterisation used to impose the constraints. However, the difference 
between the Higgs scalar masses and the squark and slepton masses must be that 
caused by the reparameterisation. Thus, the reparameterisation which makes the 
squark/slepton soft masses equal to the Higgs soft masses must be that which returns 
us to the initial state. This can be called the "universality reparameterisation" . Re- 
imposition of the high-scale boundary conditions then leaves the model unchanged. 

As the calculation converges toward a solution consistent with both sets of 
boundary conditions, we need an estimator for the difference between the Higgs and 
squark/slepton soft mass-squareds at M X - Until convergence is reached, m Hl ^ m H2 , 
rriu ^ mj etc. The choice used here is the difference between the average Higgs mass- 
squared and the average squark mass-squared. 

Unlike the /i = approach this approach allows us, for any point in the MSSM, 
to model the same point, by using the NSSB with C — jl = 0. The resulting solution 
is identical with that obtained using SOFTSUSY normally, and the high-scale value 
of \i which the model selects to fit tan (3 is the same as in the MSSM. It is possible 
to gradually turn on jl and Cq and move away from any MSSM point. The MSSM 
is a D — 2 dimensional subset of the D dimensional NSSB parameter space. 
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Table 1: The high-scale boundary condition(b.c) part of the iterative procedure for 
determining the model parameters consistent with the chosen input parameters. The state 
of the model after convergence is achieved is shown, i.e., once the model has become 
invariant under the iterative numerical procedure used. Thus, the imposition of the high- 
scale boundary conditions leaves the parameters unchanged, as shown in the table. m| 
is the average squark mass-squared, is the average Higgs mass-squared. X c is used to 
denote a particular value of the parameter named X. 

6. NSSB /3-functions 

Ref. [7] presents general one-loop beta functions for a general SUSY model with arbi- 
trary particle content and gauge group, allowing NSSB terms. The two-loop results 
are presented in [8]. Reference [7] also gives the RGEs specialised to the NSSB MSSM 
with one family and with the Yukawa couplings factored out of the parameters. For 
example, the Lagrangian "A" soft breaking term is written as mioY t H2Qt. 

An alternative presentation of the beta functions, better suited to numerical 
implementation, is the 3-family matrix form for the MSSM without NSSB used 
in [16], with the Yukawa matrices not factored out of the soft breakings. This version 
of the beta-functions has been calculated. 

A subtle question is the impact of the /z-reparameterisation on the beta functions. 
The results of [7] are in ji — reparameterisation. The MSSM results of [16] have 
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non-zero /i. We do not wish to re-calculate all of the MSSM part of the beta functions 
for the ix = reparameterisation. We wish to use part of the results of [16], modified 
according to the general results of [7], which are for fx = 0. By considering the fx- 
reparameterisation behaviour of each quantity, and representing the beta functions 
as pieces with different reparameterisation behaviours, and with dependence on each 
parameter, it is possible to obtain for the soft masses: 

% % % 

where the left hand side of each equation is the NSSB beta function for fx = 0, 
0MSSM j g MggM fj, ^ beta-function with fx set to zero, and /3 is the piece 
from the general result of [7] which has jx or C dependence. A similar equation can 
be obtained for j3 B , but it is simpler just to calculate it afresh. The beta functions 
for C and jx can be obtained from the general expression of [7]. 

The renormalisation group equations are presented in appendix 

7. Radiative electroweak symmetry breaking and other low- 
scale considerations 

The constraints given by analytical minimisation of the Higgs effective potential are 
used to obtain fx and B from tan/3 and v. This is done at the scale (Ms) where 
radiative corrections to the renormalisation group improved scalar potential have 
the smallest scale dependence, see below. 

The form of the tree-level radiative electroweak symmetry breaking (REWSB) 
constraint equation for fx: 

1 „ ,2 _ m 2 Hi +Lt 2 -(m 2 H2 + Lt 2 ) tan 2 p 

2 Z ~ tan 2 /5-l (7 ' lj 

suggests that the resulting value of fx will be correct in any reparameterisation, 
because of the way in which the reparameterisation keeps m 2 H + fi 2 constant. This 
equation finds the \x value consistent with given Higgs soft masses. However, in the 
/i = reparameterisation, numerical implementation of this equation fails, because 
the resulting fi 2 value close to convergence toward zero is sometimes negative. 

The values of the model parameters obtained after imposing the REWSB con- 
straints, if compared in the same reparameterisation, must be the same whatever 
reparameterisation we impose eq. (|7.1|) in. This is a consequence of the invariance 
of the Lagrangian, and can be verified by considering the application of a reparam- 
eterisation to the tree-level result eq. (|7.1|) . In order that the REWSB routine be 
convergent around a positive value of /i 2 , we can thus choose to impose the REWSB 
constraint after reparameterising so that jx = 0. 

The constraint equation ( |7. 1|) is obtained by minimisation of the effective po- 
tential. In our study, such a constraint equation is used to obtain fx, and a similar 
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expression to obtain B. The one-loop corrections to the effective potential are in- 
cluded, with the leading top/stop contributions from the NSSB couplings, and the 
calculation is performed at the renormalisation scale where the scale dependence of 
the 1-loop corrections is smallest, approximately ^/m t -m t - . This scale, Ms is ap- 
proximately 760 GeV for the mSUGRA parameters used here, and is intermediate 
between Mz and Mx- It is believed that the scale at which the 1-loop result has 
the slowest scale dependence, the contribution from higher loops is smallest. This is 
motivated by ref. [17], which shows that for the MSSM Higgs potential, the tree-level 
and 1-loop results are closest at this scale. However, this relationship has been shown 
not to work as well for other potentials, such as MSSM charge-and-color-breaking 
potentials, in ref. [18]. 

The physical mass matrices for those particles affected by pL and C in \x = 
reparameterisation are given in [7]. The unaffected mass matrices can be used as 
originally implemented in SOFTSUSY. 

The NSSB contributions to the mass matrices given in [7] are the leading, tree- 
level, NSSB corrections to the relationship between the physical particle mass and 
the running model parameters. SOFTSUSY includes some MSSM loop corrections 
in the mass matrices, but such corrections involving the NSSB couplings have not 
been added. 

8. Results on the particle spectrum 

As an example of the superparticle spectrum in NSSB, figures |2| and |3| show some 
of the masses on lines through NSSB parameter space passing through mSUGRA 
standard point II (defined using the sign conventions of [19]): 

m = 400 GeV; 
Mi/2 = 400 GeV; 
tan (3 = 10; 
A = 0; 

\l > 0. (8.1) 

The same point with the opposite sign of //, has been investigated and shows 
very similar results, except with Cq t— > —Cq and jl i— > —jl. 

Note that there is nothing special about the spectrum at the zero values of Cq 
or jl. The MSSM point already has broken supersymmetry. The NSSB terms break 
a few remaining equalities between SUSY breaking terms, but this has no sudden 
qualitative effect on the superparticle spectrum. 

The masses and soft breakings for nonzero Cq near mSUGRA point II show an 
extremum in fig. 0. This occurs where the B parameter is equal at Ms and Mz- 
This behaviour is a property of the RGEs for the system. 
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Figure 2: Some of the MSSM particle masses, and the B parameter at the REWSB scale, 
for nonzero values of the non-standard supersymmetry breaking parameter Cq. Cq = 
here is the mSUGRA point II: eq. (|S.l|), At the missing point, and the final point, the 
numerical method failed to converge correctly. Detail of this region can be seen in figure |^. 
The masses and parameters shown are those referred to in section |9|, or those masses which 
vary significantly across the range. Those masses shown which do not change significantly 
have been chosen to indicate this fact. 

The very low masses of the normally heavier Higgs bosons (A, H^, H°) for large 
negative values of Co, as shown in fig. |2|, and the removal of all but the lightest 
neutralinos and charginos from the lower mass part of the spectrum (fig. |]) for large 
negative values of ft, might show interestingly different phenomenology. This may 
be worth investigating with ISA JET [14] or HERWIG [20]. 

9. Parameter space constraints 

What determines the largest positive and negative values of Co and fx allowed? For 
very large negative C both the eigenvalues of the Higgs 1-loop effective potential 
at Ms become negative at (Hi) = (H 2 ) = 0. This is in contrast with the usual 
MSSM case where only one of the eigenvalues is negative - a phase transition has 
occurred. In this region, the REWSB constraint (eq. ( |7.1|) ) fails to correctly place 
the minimum of the Higgs effective potential at the required values of v and tan/3, 
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Figure 3: Some of the MSSM parameters and particle masses for nonzero values of the 
non-standard supersymmetry breaking parameter /x. /x = here is the mSUGRA point II: 
eq. Q. 

placing a saddle point there instead. See fig. ^ where the Higgs masses become small 
as we approach the transition at Co ~ —750 GeV. The shape of the surface below 
the transition is shown in fig. f|. Note that eq. assumes that the point found 

does correspond to the physical minimum, so we do not have a correct determination 
of the superspectrum in the unusual phase. 

In order to obtain the above understanding, and obtain figure |], code was written 
to obtain the full shape of the Higgs effective potential away from the minimum, 
including the dominant 1-loop corrections. 

We should next attempt to resolve the question: is there a point in the parameter 
space where this new phase yields physical values of the particle masses? For the 
parameter set shown in fig. |4] the Higgs VEV is much too large, resulting in a value of 
Mz which is too large. Perhaps, for a given pair of tan j3 and | (H) | , more than one set 
of values for the soft breakings is possible, corresponding to each of the two phases. 
Although, at tree level, eq. (|7.1|) suggests a unique value of the parameters for a given 
tan /3 and | (H) | , it must actually be solved iteratively when the loop corrections are 
taken into account, because these corrections depend on the soft SUSY breakings. It 
is therefore possible that this iterative procedure possesses more than one solution. 
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< H 1 > / GeV 

Figure 4: The Higgs effective potential surface for large values of Cq near mSUGRA point 
II. (This plot was made with Co = —800 GeV.) Note the saddle point at tan/3 = 10, (v) = 
243 GeV, i.e. at (H 2 ) ~ 240, (Hi) ~ 0. The central (red, solid) contour is -1 x 10 8 GeV 4 , 
the next contour (magenta, solid) is at —8.5 x 10 9 GeV 4 , the small elliptical contours 
(blue, dot-dashed) are at -2 x 10 10 GeV 4 , the blue, dashed contour is 10 11 GeV 4 , and the 
outermost (green, solid) contours are at 2 x 10 11 GeV 4 . 

However, after an extensive search of parameter space, nowhere was the new phase 
found to be relevant to real physics. 

There exist points in parameter space where the "universality reparameterisa- 
tion" discussed in section |4] does not exist after imposing the low-energy constraints 
and running up to the high scale: the required value of fj, 2 in the universality repa- 
rameterisation at high scale becomes negative. At such points the Higgs soft masses 
cannot, for positive fi 2 be made equal to the sfermion soft masses, while maintaining 
consistency with the electroweak symmetry breaking conditions. If this equality is 
forced, there is no positive value of fi 2 for which electroweak symmetry breaking 
occurs. This provides the limit on large negative values of jl at mSUGRA point II. 

It is for these large negative values of jl at mSUGRA point II that the recurrence 
relation approach becomes unstable. As shown in figure [5], the convergence toward 
the fixed point becomes steadily slower, before it finally becomes unstable. A naive 
belief in the recurrence relation approach would result in an underestimate of the size 
of the valid parameter space. The recurrence relation approach becomes unstable at 
jl ~ —675 GeV, while the single-ended drift method explained in appendix |B| extends 
the region which can be explored as far as jl ~ —775 GeV. Even that method cannot 
reach the physical limit, due to numerical instabilities caused by early iterations 
entering the unphysical region. This is illustrated in fig. 
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Figure 5: At mSUGRA point II, the recurrence relation solution to the SUSY RGEs 
converges more and more slowly for increasingly negative values of jl. Each separate con- 
vergence in the plot is for a different value of jt, with jt becoming more negative to the 
right. Note that for sufficiently large negative values of jl, the method fails to converge at 
all. The numbers near the lines show the value of jl used for that point, in GeV. 



For large positive jl at mSUGRA point II, the bound on jl due to the value 
of B runs to extremely large negative values at high scale, making the calculation 
unstable. This is an artificial consequence of using B rather than m\. As shown in 
fig. |3|, mj = B(fi — jt) is finite, but jj, — jl is becoming small, so B is singular. The 
calculation becomes unstable here due to the large value of B, but the upper bound 
on the value of jt comes from the lightest neutralino mass, which is becoming small 
as (/i — jl) does. 

For large positive values of Co, above around Cq ~ 1500 GeV, the drift method 
becomes non-convergent. The single-ended shooting method discussed in section ||] 
and appendix || takes us further, and is used to produce figure H, which takes us up 
to Co ~ 1675 GeV. Extrapolating, it appears that the physical limit occurs for the 
same reason as for large negative Co, i.e. the Higgs masses become small. However, 
the physical limit cannot be reached even with the single-ended shooting method, 
for the reason discussed above. With care, some useful data can be teased out of 
this region (see fig. 0), showing that the physical limit is around Co ~ 1675 GeV. 



14 



300 



250 



200 



m 



150 



/GeV 




100 



500 



Figure 6: The value of n in universality reparameterisation at Mx for non-zero /2 near 
mSUGRA point II. The dotted part of the line shows the extra reach provided by the 
single-ended drift method. 



10. Fine-tuning 

Eq. (|7.1|) shows that, with rriH 1 and m# 2 of order the mass scale of the scalar super- 
particles, a cancellation is required so that Mz can take its known value. The larger 
the typical mass of the scalar superparticles (controlled by the mSUGRA parameter 
m ) the more finely tuned this calculation must be. One possible numerical measure 
of the delicacy of this calculation is [21,22]: 



C a 



din Ml 



din \a\ 



(10.1) 



The fine-tunings with respect to each model parameter a are calculated, and the 
overall fine-tuning of the model is defined to be c = max(c a ). 

Here, there is a subtlety, again caused by the presence of the reparameterisation 
invariance: the same Lagrangian may be described by a different choice of parame- 
ters. For example, one could choose, instead of Co and /i, to reparameterise so that 
Co = 0, and use (jn^J — WTh) as the other parameter. Then, one would get different 
fine tuning values for the same Lagrangian, depending on ones choice of parameter 
set. If this situation is thought to be unsatisfactory, the solution is to take the fine 
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Figure 7: The mass of the heavy neutral Higgs particle, for non-zero Co near mSUGRA 
point II. All the data was produced using the single-ended shooting method. The dashed 
part of the line is the region near the physical limit where transients spoil the stability of 
the method. The figure shows the outcome of the method, whether or not it has formally 
converged. 



tuning with respect to the reparameterisation-invariant combinations of parameters 
appearing in the Lagrangian, for example, (p, — fi). 

On the other hand, one may view it as acceptable that the fine-tuning depends 
upon the choice of reparameterisation, given that the fine-tuning measure is sup- 
posed to be a measure of ones satisfaction with the deeper model lying beyond 
the Lagrangian being studied. A parameterisation choice reveals something of ones 
assumptions about the beyond-the-MSSM physics giving rise to the MSSM-type La- 
grangian. 

Fig. |8] shows fine-tuning contours near mSUGRA point II with respect to invari- 
ant combinations of parameters, and fig. for the parameter set chosen here. 

These results show that, in both cases, there is significant fine-tuning with respect 
to negative fb. There is no significant fine-tuning with respect to Co, and for positive 
jl the dominant fine-tuning is that with respect to the standard supersymmetry 
breakings. These effects are consistent with the behaviour explained in section || 
The absolute magnitude of the fine tuning is smaller when considered with respect 
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Figure 8: Fine tuning contours for non-standard super symmetry breaking, with fine- 
tuning differentials taken with respect to invariant combinations of parameters. The plane 
is shown parameterised in terms of Cq and ft. The outermost line shows the edge of the 
valid parameter space, according to the drift method. 
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Figure 9: Fine tuning contours for non-standard supersymmetry breaking, with fine- 
tuning differentials taken with respect to the usual mSUGRA parameters, Co and ft,. 



to the invariant combinations of parameters. 



11. Summary 

The 3-famiiy beta functions have been obtained for the MSSM with the NSSB terms, 
in a form suitable for numerical implementation, and in an economical way making 
use of the reparameterisation invariance. 
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The mSUGRA boundary conditions have been adapted for use in the context 
of the NSSB reparameterisation invariance, by choosing to impose the boundary 
conditions in the reparameterisation where they are most nearly satisfied. These 
boundary conditions allow the conventional mSUGRA (or CMSSM) model to be seen 
as a subspace of the full parameter space including NSSB terms, and the calculation 
of the superspectrum can be made in the conventional way. 

The superparticle spectrum and a fine-tuning measure have been calculated for 
non-standard SUSY breaking near mSUGRA point II, where it is found that there 
is no significant experimental signature indicating the new terms, except where they 
approach their maximum allowed values. 

The conventional MSSM methods are, however, found to be inadequate in this 
extreme region, due to the instability of the numerical method used. A more stable 
numerical method has been developed which solves this problem. 

The phase transition in the Higgs potential occurring for large negative C at 
mSUGRA point II suggests that this model potentially has a rich phase structure, 
although this appears to be of only theoretical significance, since in this phase, the 
Higgs VEVs are very large. 

In conclusion, we see that the numerical approach used to calculate the MSSM 
superspectrum can be applied to this model, but there are subtle complexities which 
must be dealt with. As supersymmetric models continue to be developed, given 
the need to calculate their properties sufficiently accurately for comparison with 
experiment, such subtleties must be treated carefully. 
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A. Renormalisation group equations 



The 1-loop three family /^-functions for the MSSM are presented below, in the repa- 
rameterisation where /i — 0, with family, colour and electroweak indices suppressed. 
The convention for the /9-functions of ref. [16] is used, i.e.: 

& = ~^-y (AT) 

where m is the renormalisation scale. 

Note that for the /9-functions of the scalar masses, only the additional contribu- 
tion from the non-standard terms is given. These should be added to the expressions 
given in [16], which is denoted as p MSSM , The p MSSM in SOFTSUSY include some 
two-loop terms. For /9-functions not given below, the NSSB couplings do not con- 
tribute, and the results of [16] may be used. 



167r 2 /5a" = 3Tr(Y d Y d ^)C u + Tr(Y e Y e ^)C u + 2C u Y u] Y u + Y U Y'^C U 
+6Y u Ti(C u Y u ^) + SY d Y d ^C u + 2C d Y d ^Y u , 

-C u {\gl + fg 2 3 ) - 4:jlY d Y d ^Y u - 2jlY u {gl + 3g 2 3 ), (A.2) 

lQ7r 2 (3 cd = 3Tr(Y u Y^)C d + 2C d Y d ^Y d + Y d Y d ^C d 

+6Y d Ti{C d Y d ^) + 3Y u Y u ^C d + 2C u Y ujf Y d + 2Y d Tr{C e Y er ), 
2 1 6 

+C\-gl - -gl) - m u Y u] Y d - 2JX\g\ + 3g 2 3 ), (A.3) 

lQ7c 2 (3 Ce = 3Tr(F M F Mt )C e + 2C e Y e ^Y e + Y e Y e ^C e 
+2F e Tr(C e r et ) + 6F e Tr(C d r et ) 

-2C e g 2 - 2liY\g\ + 3g 2 3 ) (A.4) 

16n 2 (3 mk = 1Qtt 2 (3^ m + 6Tr( Wt) - ( 6 ^ 2 2 + 2^)/i 2 , (A.5) 

lQn 2 (3 m 2 = 1Q7T 2 (3^ SSM + QTr(C d C d ^) + 2Tr(C e C ef ) - (6g 2 + 2g 2 1 )jl 2 , (A.6) 

167r 2 (3 m 2 = 16n 2 (3^2 SSM + 2[C U C U ^ + C d C d] ) - Ajl 2 (Y u Y u] + Y d Y% (A.7) 

l^ 2 (3 m 2 = 16tt 2 (3^2 SSM + 2C e C et - Ajl 2 Y e Y e \ (A.8) 

\Qn 2 (3 m 2 = \Qn 2 (3^2 SSM + 4C ut C M - 8j2 2 Y u ^Y u , (A.9) 

U R "i? 

167r 2 /3 m 2 = \Q7i 2 f3^2 SSM + AC d ^C d - 8fi 2 Y d ^Y d , (A.10) 

lQn 2 (3 m 2 = 16vr 2 /3^ 5SM + 4C et C e - 8/i 2 F et F e , (A.ll) 



e R e R 

2 



16tt 2 /3 b = -Tr(3C u U„ + 3C d ^A d + C e] ' A e ) + 2g 2 1 M 1 + 6g 2 M 2 , (A. 12) 



16tt% = jxTi{3Y u Y u ^ + 31"^ + Y e Y eX ) - g\ - 3g\. (A. 13) 
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B. Numerical methods 
B.l Introduction 

The numerical methods used in this paper may be formally described as follows, 
where for clarity we note what each general variable corresponds to in our SUSY 
case. Let the boundary conditions at scale x = Mx be 

f(p(x))=0. (B.l) 

This equation corresponds to the mSUGRA boundary conditions, e.g. m 2 q = and 
p(x) corresponds to m , M 1 / 2 , etc. Let the boundary conditions at scale z (= M z 
and/or Ms) be 

g(q(z))=0. (B.2) 

This condition corresponds to the physical values of the Yukawa couplings, {h t etc.), 
the gauge couplings, and the REWSB condition (which constrains \i and B). The 
full set of parameters for the model is {p(t),q(t)}, with p(t) constrained at t — x 
and q(t) constrained at t — z. Let qc be the q solving eq. (p.2j ) and let pc be the p 
which solves eq. ( P-l| ). The differential equations (renormalisation group equations) 
then relate the parameters at the two scales: 

q(z) = R q (q(x),p(x)) 

p(z) = R p (q(x),p(x)). (B.3) 

B.2 The shooting method 

For a given q(x), i.e. /z, B, tan /3, and #1,2,3, at Mx, therefore, g(q(z)) can be calcu- 
lated numerically through R q . The solution of the boundary value problem therefore 
reduces to the numerical solution of the algebraic equation g(R q (q(x),pc)) = for 
q(x). This is the "shooting method" of [23]. 

The shooting method cannot, however, be used for SUSY. Firstly, we have three 
scales at which boundary conditions are imposed: Mz,Ms and Mx- Secondly, the 
boundary conditions at Mz and Ms cannot be conveniently expressed in the form of 
eq. QB3). 

B.3 The drift method 

In most implementations of the shooting method for SUSY, an alternative technique, 
which can be called the "drift method", is used: Start with a guess, q 9 , for q(x), 
i.e. fi(Mx) etc, and determine p(x) from eq. ( B.l|) . Determine q(z) and p(z) from 
eq. ( B.3|) , the renormalisation group equations. Impose eq. (p.2j ) to re-determine 
q(z). (Set the Yukawa matrices and gauge couplings from experimental data and 
set fi and B from the REWSB constraints.) Use the inversion of |B.3j , i.e. run back 
up to high scale, to find new values for q(x) and p(x), q'(x) and p'(x) respectively. 
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This corresponds to a recurrence relation q'{x) = Q(q 9 ). Alternatively, the same 
procedure can be applied with a guess p 9 for p(z), i.e., the soft breaking parameters 
at Mz, and a recurrence relation for p(z). 

The required values for p, q consistent with both sets of boundary conditions 
will be a fixed point of Q. We can hopefully determine this fixed point by making 
a guess for q(x) and iterating the procedure until it is convergent, i.e. hoping the 
recurrence relation will drift to the fixed point we want. Convergence is tested for 
by requiring q and p not to change over a cycle. However, note that the fixed point 
corresponding to the physical solution is not guaranteed to be the only fixed point, 
nor is it guaranteed to be stable. 

B.4 The single ended shooting method 

When the required fixed point of Q becomes unstable, and/or the drift method drifts 
into a different fixed point of Q, a different approach can be used. The equation: 

Q( q 9) - q 9 = Q (B.4) 

can be solved for q 9 . In this work, the Newton- Raphson method was used to solve 
this equation. This has improved stability, but is not perfect, because any unrequired 
fixed points of Q will also correspond to roots of eq. (|B.4|) . The method can be slow, 
but, if an approximate solution only is desired, one can just take initial guesses for 
the less important parts of q 9 such as Y e , rather than solving for them, making the 
method faster by reducing the parameter space to be solved. 

The procedure used to generate q' also generates p' = P(q 9 ). For q 9 consistent 
with both sets of boundary conditions (and therefore a solution of eq. ( |B.4|) ), p' = p c , 
and therefore f{p') = 0. This provides a check that one has found the desired root. 

One might imagine that the boundary value problem can be reduced to the 
numerical solution of the algebraic equation P(q 9 ) — pc = 0, but this turns out to 
be unstable: formally, it is overconstrained (there being more elements in p than in 
q for SUSY), but many of the constraint equations are approximately degenerate. 
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